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4 Abstract — A tropospheric ozone Differential Absorption Lidar 

5 system, developed jointly by The University of Alabama in 

6 Huntsville and the National Aeronautics and Space Adminis- 
aqi 7 tration, is making regular observations of ozone vertical dis- 

8 tributions between 1 and 8 km with two receivers under both 

9 daytime and nighttime conditions using lasers at 285 and 291 run. 

10 This paper describes the lidar system and analysis technique 

1 1 with some measurement examples. An iterative aerosol correction 

12 procedure reduces the retrieval error arising from differential 

13 aerosol backscatter in the lower troposphere. Lidar observations 

14 with coincident ozonesonde flights demonstrate that the retrieval 

15 accuracy ranges from better than 10% below 4 km to better 

16 than 20% below 8 km with 750-m vertical resolution and 10-min 

17 temporal integration. 

18 Index Terms — Differential Absorption Lidar (DIAL), lidar, 

19 ozone, remote sensing, troposphere. 

20 I. Introduction 
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which has monitored ozone since the 1960s. The ozonesonde 36 
profiles ozone with a 100-m vertical resolution from the surface 37 
to 35-km altitude with the accuracy of 5%-10% [15], [16]. 38 
Ozonesondes are attractive because of their low up-front cost 39 
and well-characterized behavior. However, they are not suitable 40 
for making continuous measurements because of logistical con- 41 
siderations. Interesting atmospheric phenomena that vary over 42 
periods less than one day are particularly difficult to monitor 43 
using balloon ozonesondes. Satellite observations can derive 44 
total column ozone [17] and stratospheric ozone [ 1 8]— [22] 45 
and extend measurements to altitudes that are inaccessible to 46 
ozonesondes. More recently, high-quality satellite observations 47 
of tioposphenc ozone are becoming available [18], [23]— [33]. 48 
Although the satellite measurements can produce global maps 49 
of ozone, their current measurement unceitainties, along with 50 
their coarse spatial and temporal resolution, limit their ability to 51 
obseive short-teim variations in ozone. Lidars can supplement 52 
these techniques when a requirement exists loi ozone retrievals 53 
with higher temporal (from 1 min to several hours) and vertical 54 
resolution (from tens of meters to 2 km) For example, lidars 55 
of the Network for the Detection of Atmospheric Composi- 56 
don Change [34], [35] are providing long-teim observations 57 
of ozone, as well as aerosol, temperature, and water vapor. 58 
Although the up-front costs are considerably higher than for 59 
a balloon ozonesonde operation, lidars can acquire profiles 60 
continuously under both daytime and nighttime conditions. 61 
The spatial and tempoial resolution of a lidar is more than 62 
sufficient to characterize short-term ozone variations for the 63 
photochemical studies of vertical processes. 64 

Differential Absorption Lidar (DLAL) has been successfully 65 
used to measure ozone within the PBL [36], [37], the free tro- 66 
posphere [38]-[44], and the stratosphere [45]-[48] for several 67 
decades. DIAL is evolving from ground-based and airborne 68 
systems to systems that aie suitable for long-term deployment 69 
in space [49] The technique derives ozone concentrations by 70 
analyzing how rapidly the backscattered signals at two sep- 71 
arate but closely spaced wavelengths, one strongly absorbed 72 
by ozone and the other less strongly absorbed, diminish with 73 
altitude. This measurement does not require knowledge of the 74 
absolute signal intensities but, rather, only the relative change 75 
of the two signals with respect to altitude. Using electron- 76 
ically gated detection permits range-resolved measurements 77 
to a resolution as small as several meters over acquisition 78 
times of several minutes. The ozone DIAL discussed in this 79 
paper is located in the southeastern U.S. and thus provides a 80 
unique observational site within an interesting scientific area 81 
[50] to study trace-gas transport at the midlatitudes for both the 82 
polluted PBL and the free troposphere. 83 


21 ZONE IS A KEY trace-gas species within the tro- 

22 \ J posphere. On the one hand, ozone is a precursoi of the 

23 hydroxyl radical [1], which reacts with most trace species in the 

24 atmosphere. On the other hand, ozone is also a strong green- 

25 house gas influencing the climate by its radiative forcing [2]. 

26 In situ photochemistry and dynamic processes largely govern 

27 the distribution of tropospheric ozone [3] Measuring ozone 

28 variability at high spatial and temporal lesolution increases 

29 our understanding of tropospheric chemistry [4], [5], plane- 

30 tary boundary layer (PBL)-free-tropospheric exchange [6], [7], 

31 stratosphere-troposphere exchange [8]-[10]. and the impact of 

32 lightning-generated NO x on tropospheric ozone [1 1 ]— [14j 

33 Several techniques currently exist for making lange-resolved 

34 measurements of tropospheric ozone. The most common tech- 

35 nique is the balloonbome electrochemical concentration cell, 
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84 n. System Description 

85 Housed in the Regional Atmospheric Profiling Center for 

86 Discovery (RAPCD), the tropospheric ozone DIAL system is 

87 located at 34.7250° N, 86.6450° W on the campus of The 

88 University of Alabama in Huntsville (UAHuntsville) within 

89 the Huntsville city limits at an elevation of 206 m above 
AQ3 90 sea level. It is designed for measurements within the PBL 

91 and the free troposphere during both daytime and nighttime. 

92 Because of UAHuntsville’s location and occasional high tem- 

93 perature and humidity conditions, heavy aerosol pollution is 

94 sometimes present. Compared with the clean free troposphere, 

95 these aerosols require a larger dynamic range for the detection 

96 system because of their larger optical depth. Moreover, the 

97 rapid change of aerosol concentrations (e.g., due to convec- 

98 tive activity) increases the measurement uncertainty for DIAL 

99 within the PBL and lower troposphere. Judicious system-design 

100 choices and an effective aerosol correction scheme allow this 

101 system to produce high-quality ozone profiles under a variety 

102 of conditions. 


ratios at 8 km under most sky conditions. The retrieval errors 139 
due to aerosol interference are a concern in the PBL and 140 
lower troposphere. These errors are not a simple function of 141 
the wavelength separations because reducing the separation to 142 
reduce the aerosol differential backscattering will also decrease 143 
the differential ozone cross section. These errors are sensitive 144 
to the local aerosol composition, size distribution, and vertical 145 
profile. Although the aerosol interference can be lower when 146 
our online wavelength extends to the steepest part of the ozone 147 
absorption cross section, this will significantly sacrifice the 148 
maximum measurable range. Therefore, the 285-29 1-nm pair 149 
is the optimal choice to balance the maximum measurable 150 
altitude, the impact of aerosol differential backscattering, and 15 1 
the impact of solar background. 152 

B. Hardware Components 153 

Table I lists the characteristics of the RAPCD ozone DIAL 154 
system. The transmitter consists of two identical dye lasers 155 
pumped by two separate frequency -doubled Nd:YAG lasers 156 
(Fig 1). A pulse generator triggers each laser pulse with a 157 
25-ms separation between the alternate pulses. The dye lasers 158 
are software controlled to select the usei-defined wavelength. 159 
The knife-edge method [52] determines that the divergences 160 
of both UV laser beams are less than 1 mrad. A 0.75-m 161 
tuple-grating monochromator (Acton Research Corporation) 162 
indicates that the actual wavelengths of the outgoing UV lasers 163 
are 285 and 291 nm within an uncertainty of 0.1 nm. 164 

The icceiving system currently opeiates with two separate 165 
telescopes, as shown in Fig. 2. The high-altitude receiver uses 166 
a 40-cm Newtonian telescope, and the low -altitude channel 167 
employs a 10-cm Cassegrain telescope. The large telescope 168 
system routinely makes measurements from 3 to 8 km and, 169 
on occasion, measures ozone at 12 km Employing a 1.5-mrad no 
field of view (FOV), the large telescope achieves full overlap 171 
between the laser and receiver at about 3 km. Larger FOVs 172 
lower the altitude at which full overlap occurs but significantly 173 
increase solar background. The small telescope system cur- 174 
rently retrieves ozone between 1 and about 5 km with a typical 175 
FOV of 4 3 mrad The future plan is to extend the retrievals 176 
down to about 200 m with an additional altitude channel in the 177 
small telescope. The bandpass filters used to restrict the solar 178 
background for both receivers have a transmittance of 35% at 179 
285 nm and 20% at 291 nm. 180 

The detection system of the RAPCD ozone DIAL uses both 181 
photon counting (PC) and analog detection to facilitate oper- 182 
ations ovei both altitude channels. This detection combination 183 
provides the linearity of the analog signal in the strong-signal 184 
region and high sensitivity of the PC signal in the weak-signal 185 
region. An EMI 9813 QA PMT, which has been used.exten- 186 
sively for many years on a number of Goddard Space Flight 187 
Center lidar systems [53], [54], is used in the high-altitude 188 
channel, while a small Hamamatsu 7400 PMT is used in the 189 
low-altitude channel. A photodiode detects the outgoing laser 190 
pulses, which trigger both the PMT gating circuits and the Licel 19 1 
transient recorder (TR) (TR40-80, Licel Company, Germany). 192 
The Licel TR offers the advantage of increased dynamic range 193 
by providing simultaneous measurements using both analog 194 


103 A. Wavelength Selection 

104 The selection of the 285- and 29 1-nm wavelengths le- 

105 suits from the balance of the following three considei ations: 

106 1) optimizing the altitude range to make retrievals; 2) leducing 

107 the impact of the solar background during daytime operation; 

108 and 3) reducing the impact of aerosol interference upon the 

109 ozone retrieval. The DIAL wavelength selection is flexible 
no and optimized for the local ozone distribution, the absorption 

111 arising from non-ozone species, the measurement range, and 

1 12 the specific system configuration, including the output power, 

113 the telescope mirror size, and the photomultiplier s (PMTs) 

114 dynamic range. Numerous publications (e.g , [51]) discussed 

115 the optimum wavelengths for tropospheric systems. Although 

116 shorter wavelengths can provide highei measurement sensitiv- 

117 ity arising from the larger ozone differential cross section, they 

118 limit the maximum measurable range due to stronger attenua- 

119 tion of ozone absorption and Rayleigh (molecular) extinction 

120 and thus require more signal acquisition time. In addition, 

121 the shorter wavelengths require more dynamic range of the 

122 detection system and might require more altitude channels. 

123 With the current transmitter power, the online wavelength of 

124 285 nm allows us to measure ozone up to 9 km under a clear 

125 sky and 7 km under aerosol loading with a 10-min temporal 

126 resolution. Because of the significant solar background during 

127 daytime operations, we choose 291 nm as the offline wave- 

128 length. Longer wavelengths will cause a significant increase 

129 in the solar background and reduce the signal-to-background 

130 ratio. To measure both wavelength channels using the same 

131 PMT and simplify the system design, we used a bandpass filter 

132 with a central wavelength of 286.4 nm and a full width at half 

133 maximum of 11 nm whose transmittance is < 10~ 8 at wave- 

134 lengths ionger than 300 nm. For a bandpass filter, the integrated 

135 sky background over the filter bandwidth and the dark counts 

136 actually determine the background for both offline and online 

137 wavelengths. For our lidar configuration, the 285- and 29 1-nm 

138 wavelength region can provide sufficient signal-to-background 
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TABLE I 

Char xteristics of the RAPCD Ozone DIAL System 


System 

Specification 

Transmitter 

Pump laseis 

NcLYAG, 20 Hz repetition rate, 5-7 ns pulse length, 300 mJ pulse' 1 at 
1064 nm, 50 mJ pulse -1 at 532 nm 

Dye 

Rhodamine 590 and 610 


Emitted UV 4 m J pulse -1 at 285 run, divergence^ mrad 

3 mJ pulse' 1 at 291 nm, divergence '1 mrad 


Tuning range 

277 to 303 nm for the final UV output 

Receiver 

High-altitude channel Low-altitude channel 


Telescope 


Band-pass filter 
Detector 

Signal processing 


Newtonian, 40-cm diameter, f/4.5, Welch Mechanical Designs 
1 .5-mrad FOV Cassegrain, 1 0-cm diameter, 

<72.3, 4.3-mrad FOV 

Center wavelength at 286.4 nm v ith a 1 l-nm 1 WHM. Transmittance is 
35% at 285 nm and 20% at 291 nm 

Electron Tubes 9813QA, about 28% Hamamatsu R7400U-03, about 
quantum efficiency 20% quantum efficiency 

LICEL Transient Recorder (TR40-80) 250-MHz maximum 
photoncounting rate, 12-bit and 40-MI Iz analog-to-digital converter, 
25-ns range resolution 



Fig, 1 . Transmitter diagram. 


195 detection and PC. The Licel TR’s highest temporal resolution 

196 is 25 ns, corresponding to a fundamental range resolution of 

197 3.75 m. It is necessary to gate the high-altitude channel off 

198 for the first 10-15 /is and the low-altitude channel for the first 

199 1 fus to maintain the PMTs linearity and minimize the impact 

200 of signal-induced bias (STB) on the background count rate. 

201 III. Data Processing 

202 A. Raw Data Processing 


signals. Foi PC at high counting rates, a second pulse arriving 208 
at the discriminator before it has recovered from the previous 209 
pulse will not be counted — a period known as dead time [55]. 210 
Experiments with a function-generator-driven LED determine 21 1 
this time to be 10 ns for the high-altitude channel and 4 ns 212 
for the low-altitude channel. Our results show that the system 213 
dead time obeys a nonparalyzable model following a simple 214 
relationship, as in (1)[56], between the true count rate Ct and 215 
measured count rates Cm, allowing the impact of dead time Td 216 
on the data to be removed 217 


203 Several operations, designed to improve the measurement 

204 precision, occur before the ozone retrieval. First, average the 

205 signal returns over 10 min and 150 m. The temporal resolution 

206 of the retrieval can be varied depending on the signal-to-noise 

207 ratio (SNR). Second, apply a dead-time correction to the PC 


Ct — 


Cm 

1 — CuTd 


( 1 ) 


Third, remove the signal background. The last 10 /us (400 218 
fundamental bins) of signals ranging up to 30.72 km (far-range 219 
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Fig. 2. Diagram of the receivers and detectors. 

220 limit), which are considered to be the background region where 

221 no laser signal returns are expected, are averaged to give an 

222 approximate background. Fourth, merge the paiallel analog and 

223 PC signals into a single profile [57] after removing the offset 

224 between the analog and PC signals [58]. We found this offset to 

225 be about 250 ns for our system by carefully comparing returns 

226 derived with clouds on both the analog and PC channels. The 

227 merged region requires that the ratio of PC to analog signals is 

228 constant. Ratios that are not constant suggest eithei an incorrect 

229 background subtraction or a wrong dead-time correction The 

230 merging threshold of the PC signal is typically 20 MHz for 

231 the Hamamatsu PMT employed in our low-altitude channel 

232 and 20-30 MHz for the EMI PMT used on the high-altitude 

233 channel. Because DIAL retrievals depend on the quality of 

234 both 285- and 291-nm signals, we combine the PC and analog 

235 signals approximately at the same altitude for both lasers to 

236 minimize the retrieval error due to the merging. Examples of 

237 the ratio of PC to analog signals and their merged region for 

238 the 285-nm signal are shown in Fig. 3. The merging threshold 

239 is 20 MHz for both altitude channels. The fifth step involves 

240 smoothing the signals to reduce random noise. Our configura- 

241 tion currently employs a five-point (o x 150 — 750 m) running 

242 average applied to returns from all altitudes; smoothing reduces 

243 the effective vertical resolution to 750 m. 

244 After initial processing, an exponential-fit correction re- 

245 moves SIB from the signal returns. This bias, caused by intense 

246 light returns from the near range (also called signal-induced 

247 noise), appears as a slowly decaying noise source superimposed 

248 on the normal returns. The causes of the SIB are related to the 
AQ4 249 regenerative effects such as dynode glow, after-pulsing effect, 

250 glass-charging effect, shielding effect, and helium penetration 

251 [59]. SIB varies widely with different PMTs. For our case, the 


SIB of the EMI 9813 is larger than that tor the Hamamatsu 252 
7400. SIB can persist for several hundreds of microseconds and 253 
can exeit a strong influence on data at the lidar’s upper range 254 
wheie both signal and noise counts become compaiable. With 255 
uncoirected SIB, the raw signal falls off more slowly at higher 256 
altitudes, resulting in lower retrieved ozone values. SIB usually 257 
has more influence on the shorter wavelength channel, which 258 
falls off more rapidly with altitude Unless a mechanical shutter 259 
physically blocks the optical path to the PMT to eliminate SIB, 260 
a model must charactenze its behavior. Cairo et al. [60] and 261 
Zhao [61] have successfully used a double-exponential function 262 
for this purpose However, this collection increases measure- 263 
ment uncertainties because both the scaling and exponential 264 
lifetimes aie difficult to determine without additional indepen- 265 
dent measuiements. A more practical technique is to employ 266 
a single-exponential fit to the residual background [42], [43], 267 
[62]. For the high-altitude channel, the function’s coefficients 268 
aie automatically determined using a single-exponential least 269 
squares fit to data acquired approximately from 100 to 160 fis 270 
after data acquisition starts where the SIB becomes dominant. 271 
The start and length of the exponential fit vary with different 272 
channels (either wavelength channels or altitude channels), 273 
atmospheric structures, and lidar configurations because these 274 
parameters affect the intensity of the detected signal. For our 275 
low-altitude channel, the SIB is weaker than that of the high- 276 
altitude channel because of the different PMT and weaker 277 
signal. However, it is difficult to automatically determine the 278 
fitting function for the low-altitude channel signal using the 279 
least squares fitting method, particularly for the 285-nm sig- 280 
nal, because the far-range signal after background correction 281 
is not completely characterized by an exponential function 282 
[Fig. 3(b)]. It is useful to optimize the exponential fitting 283 
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(a) (b) 



Fig. 3. Examples of signal merging and SIB correction tor tht 285-nm signa‘ The 10-min averaged data occurred at 13 00 local time on October 18, 2008. 
(a) Normalized ratio of PC to analog after background and dead-time corrections for the low-altitude channel signal, (b) Comparison ot the non-SIB-corrected 
signal, the SIB-corrected signal, and the model as well as the SIB filling tunction, for the low-altitude signal The model uses the coincident ozonesonde 
measurement assuming no aerosol. The SIB fitting function (exp(— 1.3 — alt • 2 • 10 _1 )) was empirically derived using previously retrieved data and coincident 
ozonesonde measurements, (c) Same as (a) but for the high-altitude channel, (d) Same as (b) but for the high-altmide channel. The coefficients of the SIB fitting 
function result from an empirical single-exponential least squares fit to the signal acquired from 100 to 160 »*s aftei data acquisition starts. 


284 function for the low-altitude channel using previous retrieval 

285 data and compare the data with coincident ozonesonde profiles. 

286 The slope of the logarithm of the SEB fitting function remains 

287 for a particular configuration (i.e., outgoing power) and could 

288 slightly change for different configurations. Those retrievals 

289 corrected using the empirically derived exponential function 

290 agree with ozonesonde profiles up to 5 km within 5% bias. 

291 Fig. 3 shows the typical effect of the SIB correction and the 

292 comparison of the fully corrected signal and the model for the 

293 285-nm signal. The model simulation employs the coincident 

294 ozonesonde measurement assuming no aerosol. 


295 B. DIAL Retrieval 

2% Excellent discussions concerning the DIAL technique occur 

297 in the publications by Measures [63], Kovalev and Eichinger 

298 [64], and Browell etal [39]. The average ozone number density 

299 ri(r+Av/2) between range r and r -j- Ar can be expressed as 

300 the summation of the signal term n \ r+ ^ v / 2 y differential 


backscattering term An* r _ | _ Ar/2) , and the differential extinction 301 
term A^ r+At/2) 302 

n (r+ Ar/2) ~ n (r+Ar/2) ^ n (r+Ar/2) + ^ n {r ' Ar/2)' (2) 


One can write the discrete forms of the three terms at the right 303 
side as follows: 304 


n (r-|-A' /2) 
Ar/2) 
^ n (r+Ar/2) 


1 i n f ^Qn(0~^off(r+Ar) \ 

2ArAaos \Pofi(r)Pon(r+Ar) J 1 * 

1 1 ( flm(r) Aff(r+Ar) \ .. 

2ArAao3 V3oh ( ,)^on(, + A r)J V 

Aoq 3 **off(r+Ar 2)) 


where the subscripts “on” and “off’ represent the online 305 
(285 nm) and offline (29 1 nm) wavelengths, respectively, P is 306 
the detected photon counts, ,6 is the total backscatter coefficient, 307 
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(a) 


0.0 0.5 1.0 1.5 

0 g (10 12 molec cm” 3 ) 

(b) 


2.0 


Fig. 4. Example of a joined ozone retrieval for the lidar data in Fig. 3. (a) Separate retrieval*, of the two altitude channels. The error bars represent the one-sigma 
statistical uncertainties. The gray envelope represents ±10fr uncertainty of the coincident ozonesonde profile, (b) Joined DIAL retrieval from the two altitude 
channels and its combined one-sigma statistical uncertainty. 


308 a is the total extinction coefficient excluding ozone, and Acro 3 

309 is the differential ozone absorption cross section. P, and 

310 a are dependent on r and the wavelength. Strictly speaking, 

31 1 A<to 3 is r dependent as well because it is a function of tempera- 

312 ture, which varies with r. By ignoring the differential scattering 

313 and extinction from non-ozone species, the DIAL equation re- 

314 duces to only n s . A n b arises from aerosol differential backscat- 

315 tering. A n c consists of differential Rayleigh extinction, aeiosol 

316 extinction, and non-ozone gaseous absorption, including O 2 . 

317 SO 2 , and NO 2 . Measurements from a meteoiological sounding 

318 can usually correct Rayleigh effects. We correct the aerosol 

319 effects when they are significantly enough, particulaily in PBL. 

320 The aerosol correction discussion appears in Section ffl-D. 

321 C. Joining Retrievals From Two Adjacent Altitude Channels 

322 Final retrievals result from joining the data from two altitude 

323 channels with a weighted average. We choose to join the final 

324 ozone retrievals instead of the raw signals because the SNRs of 

325 the two altitude channels at the joining altitude are significantly 

326 different. If the retrievals derived from two different channels 

327 are statistically independent, the best estimate of these measure- 

328 ments is the two-channel weighted average [65] 

2 2 

73-best ~ ^ ^ /S> (6) 

i=l s=l 

329 where ti, is the ozone retrieval of channel i and the weight tr* 

330 is the inverse square of the corresponding statistical uncertainty 

331 ( £u , which will be discussed in Section V) 


Typically, the low- and high-altitude channels join between 
3.3 and 4 4 km. Fig 4 shows an example of a joined ozone 
profile, as well as the combined one-sigma statistical uncer- 
tainties. 

D. Aerosol Correction 

In a polluted area, aerosols can be a dominant eiror source 
in the lower troposphere. Based on (4) and (5), the vertical 
gradient of aerosol backscattering determines An 6 , and the 
magnitude of the differential aerosol extinction coefficient de- 
termines An e . The largest aerosol correction usually occurs in 
an inhomogeneous aerosol layei (i.e., the top of the PBL). One 
can solve for the ozone and aerosol profiles simultaneously with 
only two wavelengths by assuming appropriate Angstrom expo- 
nents and constant hdai latios [66], [67]. If a third wavelength 
is available and is close to the DIAL wavelength pair, one can 
use the dual-DIAL technique [68], [69] to reduce the error due 
to aeiosol When the third wavelength is far from the DIAL 
wavelength pair, one can use the method suggested by Browell 
et al [39] to correct the aerosol interference. Without the third 
wavelength, we employ an iterative procedure to retrieve ozone 
and conecl aerosol effects. To illustrate this method, start with 
the equation for ozone number density using only the 291-nm 
signal [63] 

^(fH-Ar/2) 

1 

2(703 Ar 


Wi = l/£ii- 

(7) 

x iln ( ^ In 

+ 

to 

1 

332 The uncertainty of nbcst is 


/\ s in 1 j in 

{ \P{y+Ar)J 

{^ + A r) +0fr + Ar ) )/(r + *r)l_ 


^lbrst 



( 8 ) 


-2 


(afr+Ar/2) 


1 

' a (r+Ar/2) 



(9) 


333 

334 

335 

336 


337 

338 

339 

340 

341 

342 

343 

344 

345 

346 

347 

348 

349 

350 

351 

352 

353 

354 

355 
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356 where gq 3 is the ozone absorption cross section, 0^ and 

357 0 are the molecular and aerosol backscatter coefficients at 

358 range r, respectively, and ct ( A / +Ar/2) and a£ +Ar/ .y represent 

359 the average molecular and aerosol extinction coefficients, re- 

360 spectively, between range r and r + Ar. The subscript 291 

36 1 is omitted for brevity because all backscatter and extinction 

362 parameters correspond to 29 1 nm. Solving for i3fo , (9) becomes 


be different theoretically. A;i£ r+Ar ^ and An^ r+Ar) can be 391 
approximated as [39] 392 




(4— 7?)AA / g( r) #(,-+Ar) 


A n\ 


2ArAcro3A 0 ff \1 + £( r ) l + ^(r- fAr) 

AA 




(14) 


(r+Ar) ' 


A of f 


(^fr+Ar/2) +4«fr'+Ar/2)) 05) 


Pfr) - PX P 




2r i(r+Ar/‘2)^03Ar 


2 


(a 


Asf 

(r-l-Ar/2) 


-h c * : ( r + Ar /2) 



where 77 is the Angstrom exponent, AA is the wavelength 393 
separation, and B ^ is the aerosol-to-molecular backscatter 394 
ratio at the offline wavelength defined as 395 

B(r) — P(r)/P(r)- 06) 


r2 (^(H-Ar) + ^(r+Ar)) 

( 7 * + Ar) 2 



( 10 ) 


363 Assuming that the lidar ratio (aerosol extinction-to-backscatter 

364 ratio), i.e., S — a A /0 A , is known for the 291-nm signal and 

365 further assuming that 


a 


A 

(r+Ar. 2) 


a; 


(r+Ar) 


- set 


’(r+Ar) 


01 ) 


366 ( 10 ) only contains the following two unknown variables the 

367 aerosol backscatter coefficient a^. +A t .) and the ozone number 

368 density n( r+ ^ r/ - 2 )- Molecular backscatter and extinction can be 

369 computed from nearby radiosonde data or from climatology. 

370 For the first iteration step, n( r+ ^ r / 2 ) can be computed horn 

371 (3) and inserted into (10). By assuming a start value 4(rcf) at a 

372 reference range and a constant S with range ,3^ can be solved 

373 by (10). Then, the first profile is substituted back into ( 10 ) 

374 to compute the second estimate by using a more accuiate form 

375 for a$ +Ar /2) as 

< + Ar/ 2) - 5 ('3(t +Ar) + C) /2 02) 


376 where d^ represents the value from the first estimate. With 

377 several iterations of (10) and (12) (we name this iteration the 

378 “aerosol iteration”), we can get a stable solution for 0 A r y which 

379 does not change significantly from one iteration step to the next. 

380 The aerosol iteration stop criterion is defined as < £ A in . 

381 is the relative total difference of the backscatter coefficients 

382 between two adjacent iteration steps and is defined as 

i ^ re ^ 

(#> 7” X) < 13 > 

£ fr,„ — 

r=r, 

383 where / represents the iteration step, is the starting range 

384 of the lidar retrieval, and ^ are the backscatter coefficients 

385 at range r and iteration step L £^ in is typically 0.01 for our 

386 aerosol retrievals. Aside from the number of iterations 

387 required for a stable solution is also related to the range res- 

388 olution of the signal. For simplicity, we assume that the power- 

389 law dependences with wavelength for the aerosol extinction 

390 and backscatter coefficients are the same although they can 


The estimate for the aerosol -corrected ozone number density 396 
profile is then substituted into ( 10 ) to calculate an updated 397 
aerosol backscatter profile, which, in turn, is used to compute 398 
an updated aerosol-corrected ozone profile. This iteration is 399 
named “ozone iteration” to be distinct with the coupled aerosol 400 
iteration process A similar iteration stop criterion, f ^ 3 < £^f n , 401 
as the aerosol iteration, can be defined for the ozone iteration 402 
by replacing the backscatter coefficient in (13) with the ozone 403 
number density. Typically, only two ozone iterations are re- 404 
quired when f °f n is set equal to 0.00 1 . 405 

The lidai ratio ( 5 ) exhibits a wide range of variation with 406 
different aerosol iefractive indexes, size distributions, and hu- 407 
nudity [70] The S measurements have been made most fre- 408 
quently at 308 [71] and 355 nm [72], [73]. The S for our DIAL 409 
wavelengths was assumed to be 60 si " 1 [74] constant over the 410 
measurement range for typical urban aerosols The Angstrom 41 1 
exponent (rj) is often seen as an indicator of aerosol particle 412 
size Values greater than two correspond to small smoke parti- 413 
cles, and values smaller than one correspond to large particles 414 
like sea salt [75], [76]. Most of the reported v\ for tropospheric 415 
aerosol are measured at wavelengths longer than 300 nm with 416 
a variation from zero to two [77], [78]. Considering that rj 417 
could be relatively small when it is applied in the UV region, 418 
we assume that 77 =■ 0.5 at our DIAL, wavelengths for urban 419 
aerosols [79]. 420 

Simulations were conducted to investigate the aerosol cor- 421 
rection in the DIAL retrieval under an extremely large aerosol 422 
giadient condition by assuming the aerosol, molecular, and 423 
ozone extinction profiles at 291 nm shown in Fig. 5. The 424 
hypothetical aerosol profile includes the following three basic 425 
regimes homogeneous, increasing, and decreasing extinction. 426 
The aeiosol extinction coefficients are set equal to 1CT 5 m -1 427 
below 1 .2 km and above 3 km to represent a background value. 428 
The resulting steep gradient between the low background and 429 
high aerosol value provides an extreme test for the aerosol cor- 430 
rection algorithm. The molecular extinction profile is derived 431 
from the 1976 U.S. Standard Atmosphere [80]. The assumed 432 
ozone extinction profile is constant with altitude and is based on 433 
a number density of 1.5 x 10 12 molec • cm ' 3 and an absorption 434 
cross section of 1.24 x 10 _1S cm 2 • molec " 1 at 291 nm [81]. 435 

Fig. 6 shows the comparison of the ozone retrieval both 436 
with and without aerosol correction, as well as the calculated 437 
aerosol profile, at 291 nm. This example calculation assumes 438 
that 77 0.5 and 5 60 sr -1 are known exactly, and there 439 
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Fig. 5. Aerosol, molecular, and ozone extinction coefficient profiles at 291 tim ^*8- ^ zone retrieval using different Angstrdm exponents (tj — 0, 0.5. 

for a model calculation of extreme aerosol effects. and *** lidar ratl0S = 40 ■ 60 > and 80 > in the 8610501 correction. 


Aerosol ext. coeff. (10~ 3 m *) 


c 

1 

2 3 

4 

t 

• • ■ ■ 1 

Aerosol model 



Calculated aerosol A A 



03 model — 



03 w/o aerosol corr. 

3 


03 with aerosol corT; 






A 1 



! A 

a 

: 

! A 

M 0 

r 

1 A 


r 

i 4 “ 

4-> 


A 

<3 






i 


L- ^ j 


i 



i 



i 

1 

0 

L . . . ! 

1 

1 ....... * 

c 

) 1 

2 3 


03 (10 lz molec cm 3 ) 

Fig. 6. Comparison of the simulated ozone retrieval without aerosol correc- 
tion against that with aerosol correction using the iterative procedure. The 
Angstrdm exponent (r?) and lidar ratio (5) were assumed to be exactly 
known at 0.5 and 60 sr~L respectively, for the aerosol correction. The aerosol 
correction dramatically improves the ozone retrieval. 

440 is no signal measurement error. With a range resolution of 

441 150 m, two ozone iterations produce the final aerosol-corrected 

442 ozone retrieval by setting = 0.001. In the process of cal- 

443 culating the aerosol profile, aerosol iterations produce a stable 

444 aerosol solution by setting — 0.01, which is approximately 

445 identical to the model aerosol profile. The aerosol correction 

446 procedure reduces the retrieval errors from ±50% to about 

447 ±5%. The residual errors are due to the numerical integration 

448 and the approximation of (14) and (15). The quality of this 

449 iterative procedure depends on the choice of S and According 

450 to (10), (14) and (15), S affects the aerosol profile retrieval, 

451 while rj affects only the final ozone correction. 

452 Fig. 7 shows the sensitivity test for S and rj in the aerosol 

453 correction assuming that 5 — 60 and — 0.5 are the correct 


values. Inaccuiate estimates of S or 77 can yield retrieval errors 454 
up to about 20 %. Larger 77 will overestimate An c , which 455 
produces less ozone, and vice versa. 77 has a smaller impact 456 
on Ar* 6 relative to A n e due to the 4 — tj factor. The impact 457 
of S is larger in the inhomogeneous aerosol layer than in the 458 
homogeneous layer. The peak error is larger for underestimated 459 


S relative to overestimated S [82]. 460 

We summarize the iterative procedure as follows. 461 

1) Calculate the first estimate of the ozone concentration 462 

from (3). 463 

2) Substitute the first estimated ozone into ( 10) to derive the 464 

aerosol backscatter profile foi the offline wavelength, and 465 
iterate to obtain a stable solution with (12). 466 

3) Calculate the differentral aerosol backscatter and extinc- 467 

tion corrections to obtain a second estimate of ozone 468 
using (14) and (15). 469 

4) With the second ozone estimate, go back to step 2. 470 

IV. Measurements 471 


Fig. 8 shows an ozone DIAL retrieval for 15 consecutive 472 
hours from 12:56 local tune, August 9, to 03 : 56, August 10, 473 
2008, with 10-min temporal integration (12000 shots) and 474 
750-m veitical range resolution using the data processing de- 475 
scribed in the previous section. The aerosol correction was 476 
made only at altitudes between 1 and 4 km using the data 477 
from the low-altitude channel because of the negligible aerosol 478 
effects above 4 km. The aerosol time-height curtain [Fig. 8(a)] 479 
exhibits moderate aerosol activity below 2 km with expected 480 
diurnal PBL variation and shorter timescale fluctuations due to 481 
PBL processes. The maximum aerosol correction in Fig. 8(b) 482 
corresponds to an ozone adjustment of 3-4 ppbv and occurs 483 
between 1.5 and 2.5 km for the largest vertical backscatter 484 
gradient. The retrievals for the two altitude channels overlap 485 
between 3.3 and 4.4 km to produce the final ozone profiles 486 
[Fig. 8(c)] that agree well with the colocated ozonesonde (EN- 487 
SCI model 2Z with unbuffered 2% cathode solution) launched 488 
at 13:49 local time. The time-height curtain of ozone’s 489 
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Fig. 8. Ozone DIAL retrievals made on August 9-1 0, 2008. (a) Calculated aeio>ol extinction coefficient at 291 nm. The feature at 2 km, 14:00 is a cloud, 
(b) Aerosol correction for ozone DIAL retrieval, (c) Ozone DL\L retrieval after aciosol correction The retrieval was made with a 750-m vertical range resolution 
and a 10-min temporal resolution. The colocated ozonesonde marked by a triangle launched ai 13: 49 local time. 


490 evolution shows a very interesting structure of multiple ozone 

491 layers in the lower atmosphere that varies with time. One can 

492 see the buildup and decay of various layers throughout this 

493 12-h period. The high-frequency variation in the high-altitude 

494 channel (> 6 km) results partly from lower SNR and higher 

495 uncertainty of the SIB correction, both of which increase with 

496 altitude. Fig. 9 shows the mean ozone profile and one-sigma 

497 standard deviation for the 10-min vertical profiles between 

498 12:56 and 15:06 local time in Fig. 8, as well as the coinci- 

499 dent ozonesonde measurement. The high-altitude channel has a 

500 standard deviation increasing with altitude due to the statistical 

501 error distribution. Its standard deviation is less than 13 ppbv 


below 8 km and increases to about 45 ppbv at 8.5 km where the 502 
285-nm lasei does not have sufficient SNR for ozone retrieval; 503 aqs 
therefore, we terminate the retrievals at 8 km in Fig. 8. The stan- 504 
dard deviation of the low-altitude channel retrievals is less than 505 
5 ppb\ below 4 km and reaches 8 ppbv at 5 km due to lower 506 
SNR. The standard deviation at 2 km is a little larger than the 507 
surrounding altitudes possibly because of larger ozone fluctu- 508 
ations or larger uncertainties of the aerosol correction in the 509 
ozone retrieval at the PBL top. The two altitude channels have 510 
consistent mean retrievals in the overlap region with discrepan- 51 1 
cies less than 5 ppbv and similar standard deviations at 3.3 km 512 
which most likely reflect the true ozone short-term variations 513 
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Fig. 9. Mean ozone mixing ratio and one-sigma standard deviation for the 
10-min vertical profiles between 12 :56 and 15:06 local time in Fig. 8 . The 
colocated ozonesonde was launched at 13:49 local time. The large error bar 
(~45%) at 8.5 km identifies the high-altitude limit of the retrievals (8 km). 

514 above the PBL as shown in Fig. 8 . The mean retrievals agree 

515 with the ozonesonde measurement within about 10 ppbv and 

516 have higher biases at the upper altitudes. 


517 


V. Error Analysis 


518 We divide the error budget of the DIAL retrieval into the 

519 following four categories: l) statistical uncertainties e\ arising 

520 from signal and background noise fluctuations, 2) errois ti 

521 associated with differential backscatter and extinction of non- 

522 ozone gases (O2, SO2, NO2, etc.) and aeiosols; 3) errors £3 

523 due to uncertainties in the ozone absorption cross section; and 

524 4) errors £4 related to instrumentation and electronics. £\ is a 

525 random error; e 2 , £ 3 , and e 4 are systematic errors. e\ can be 

526 written as 141] 


£1 = 


2nArA<7o 3 Vf? (SNR, X ) 2 ' 


(17) 


527 With the assumption of a Poisson distribution go\erning PC, 

528 the SNR at wavelength A and range registration j becomes 


SNR jiX = 


PjA 

(Pj.. x + Pb + Pd ) 1 ' 2 


(18) 


529 where is the solar background counts and Pd is the dark 

530 counts. It is straightforward to show that ei is proportional 

531 to (A^ATAPl) - 1 / 2 , where N represents the total number of 

532 shots, A is the unobscured area of the telescope’s primary 

533 mirror, and Pl is the number of emitted laser photons. Ar 

534 must be chosen large enough to produce an acceptably small 

535 error. Fig. 10 shows the estimated statistical errors for the 

536 high- and low-altitude channels for a 10-min integration and a 

537 750-m range resolution. e\ is typically less than 10 r /c below 

538 4 km for our low-altitude channel and could be 20% at 5 km. 

539 This altitude performance gives us sufficient overlap for the 

540 two altitude channels under most atmospheric conditions. In 

541 the high-altitude channel, e\ exceeds 25% of the retrieval ozone 

542 near 8 ± 1 km, where we terminate the retrieval. 



Fig. 10. Estimated statistical errors for the high- and low-altitude channels 
using 10-min integration and 750-m range resolution. The nighttime and 
daytime statistical errors are modeled by using the annually averaged local 
ozonesonde profile the 1976 U.S. Standard Atmosphere, an urban aerosol 
model [83], and the lidar parameters in Table 1. The ozone profile used for 
summer daytime errors is assumed 20 % higher than the annual average. 

£2 includes the interference from 0 2 , S0 2 , N0 2 , air mole- 543 
cules, and aerosols Table II summarizes the potential errors 544 
in the DIAL retneval for 285- and 291-nm wavelengths due 545 
to non-ozone absorption gases [84]— [88] The calculation of 546 
the oxygen dimer (O 2 -O 2 ) interference includes some un- 547 
certainties due to the absorption cross-sectional measurement 548 
The O 2 -O 2 absorption theory has not been entirely established 549 
[89]. Local S0 2 and N0 2 profiling data are not available. How- 550 
ever, the estimated error due to either S0 2 or N0 2 using the 551 
latest ground observation is less than 1 %. The impact caused by 552 
differential Rayleigh extinction results in an inaccuracy of less 553 
than 1% using balloon ozonesonde tetnevals of atmospheric 554 
density or by employing climatological models. 555 

The main concern comes from the aerosol interference, 556 
which depends on both the wavelengths and wavelength sep- 557 
aration. Although the aeiosol optical properties could be re- 558 
trieved from a thud wavelength, the differential effect for a 559 
DIAL wavelength pair still has some uncertainty due to the 560 
assumption for lidai idtio and Angstrom exponent. Within the 561 
PBL, where the statistical errors are small, differential aerosol 562 
backscattenug and extinction dominate the error sources [39], 563 
[41]. [43] However, it is reasonable to believe that the error 564 
due to aerosol interference is smaller than 20% after the aerosol 565 
correction, as shown in Section III-D. 566 

The uncertainty in the Bass-Paur ozone cross sections is 567 
believed to be less than 2% [81], [84], [89]. £3 will be less than 568 
3% after considering the temperature dependence. 569 

64 could be caused by a misalignment of the lasers with 570 
the telescope FOV, imperfect dead time, or SIB correction. 571 
Dead time distorts the near-range signal, and SIB distorts the 572 
far-range signal. Because the dead-time behavior is reliably 573 
characterized, the error caused by SfB usually is larger than 574 
the dead-time error. These errors related to the signal non- 575 
linearity can be experimentally diagnosed by a function- 576 
generator-driven LED laser simulator [ 90], [91]. For the 10-min 577 
integration data, £4 is estimated to be < 5% at 1-4 km for our 578 
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TABLE II 

dial Retrieval Errors Due to Non-Ozone Absorption Gases 


Gases 

Aa, differential 

References 

Mixing ratio 

References for 

O 3 retrieval 


absorption cross-section 
(cm 2 molec" 1 ) for 285 and 

for Act 

(ppbv) 

mixing ratio 

error (%) 


291 nm 





o 3 

1.15X10-' 8 

Bass and 
Paur 1981 

60 





[84] 




0/ 

4.5xl<r 27 

Fally et al. 
2000 [85] 

2.1x10 s 


1.5% 

S0 2 

-4.8x1 O' 20 

Rufus et al. 

13 b 

NREM 2006 

-0.9% 



2003 [86] 


[88] 


no 2 

-2.25x1 O' 20 

Bogumil et 

18 c 

NREM 2006 

-0.6% 



al. 2003 [87] 


[88] 


Total 





±1.5% 

a due to O 2 -O 2 


b maximum 24-hr average in 1994. Latest local monitoung data available. 
c Annual arithmetic average in 1993. Latest local monitoring data available. 


TABLE III 

Summary of the Errors in RAPCD Ozone DIAL Measurements 


Errors 

Low-altitude channel 
(1-4 km) 

High-altitude channel 
(3-8 km) 

1 . e j , statistical error 

<10% 

"25% 

2. e 2 , interference by non-ozone species 
Aerosol 

<20% 

<5% 

Non-ozone absorption gases 


<1.5% 

Rayleigh 

-1% using local radiosonde profile 

3. £ 3 , due to uncertainty in A Oq~ 


0% 

4. e 4 , due to SIB and dead-time 

<5% 

<10% 

Total RMS error 

<23% 

-"28% 


* The errors are estimated by assuming a 60 ppbv constant ozone mixing ratio ui the tropospheie 
for data with a 750-m veitical resolution and 10-min integration. 


579 low-altitude channel and < 10% for our high-altitude channel 

580 below 8 km based on our LED test results and the analysis of 

581 our previous data such as Figs. 8 and 9. A summary of the errors 

582 in the DIAL measurements is shown in Table III for a constant 

583 tropospheric ozone of 60 ppbv, 750-m vertical resolution, and 

584 10-min integration. 

585 Fig. 11 shows a comparison of 12 lidar retrievals and their 

586 single coincident ozonesonde measurement between 13 : 00 and 

587 14 : 00 local time except for the first profile on August 17, 2008 

588 (upper right panel), which was taken at 08 : 00. The aerosol 

589 correction was made at altitudes between 1 and 4 km by setting 

590 the reference altitude at ~6 km and = 1.67 x 10 7 m 1 • 

591 sr _1 [83]. Fig. 12 shows the mean percentage differences and 

592 their standard errors of the mean for all those retrievals. The li- 

593 dar retrievals of the low-altitude channel agree with ozonesonde 

594 measurements within 10% from 1 to 4 km. The relatively 

595 high errors at about 2 km possibly relate to residual aerosol 

596 correction errors around PBL height. The lidar retrievals from 

597 the high-altitude channel agree with ozonesonde to within 20% 

598 below 8 km. The statistical error and the uncertainty associated 


with the SIB collection result in larger errors for the high- 599 
altitude channel above 6 km. 600 

VI Conclusion and Future Plans 601 

The RAPCD ozone DIAL system measures tropospheric 602 
ozone piofiles during both daytime and nighttime using the 603 
285-/29 1-nm wavelength pair. The low-altitude receiving chan- 604 
nel makes ozone measurements at altitudes between 1 and 5 km 605 
using a 10-cm telescope and Hamamatsu R7400U PMTs. The 606 
high-altitude channel measures ozone between 3 and about 607 
8 km using a 40-cm telescope and EMI 9813 PMTs. Model 608 
calculations demonstrate that the iterative aerosol correction 609 
procedure significantly reduces the retrieval error arising from 610 
differential aerosol backscatter in the lower troposphere where 611 
the quality of the aerosol correction depends on the accuracy of 612 
the a priori lidar ratio and Angstrom exponent. A comparison 613 
of the lidar retrievals and coincident ozonesonde measurements 614 
suggests that retrieval accuracy ranges from better than 10% 615 
after the application of an aerosol correction below 4 km to 616 
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O-(10 12 molec cm" 3 ) 


Fig. 1 1 . Comparison of the (solid) low- and (dashed) high-altitude-channel aeroso^corrected retiievals with the (dotted) coincident o/ontsonde measurements. 



Mean percentage difference (%) 


Fig. 12. (Dark) Mean percentage differences, lidar sonde/sonde, and (gray) 
their estimated one-sigma standard error of the mean for the data in Fig. 1 1 . 

617 better than 20% for altitudes below 8 km with 750-m vertical 

618 resolution and 10-min integration. Error sources include sta- 

619 tistical uncertainty, differential scattering and absorption from 

620 non-ozone species, uncertainty in ozone absorption cross sec- 

621 tion, and imperfection of the dead-time and SIB corrections. 

622 The uncertainty in the SIB correction and the statistical errors 

623 dominaie the error sources in the free troposphere and could be 

624 reduced by increasing the integration time or reducing the range 

625 resolution. 

626 Future improvements will overcome two major limitations 

627 of the current system by doing the following: 1) extending 

628 observations into the upper troposphere by replacing the current 

629 transmi.ters with more powerful ones and shifting the current 


wavelengths to longer ones to make higher-altitude nighttime 630 
measurements and 2) minimizing aerosol interference in the 631 
lowei troposphere by adding a third wavelength (dual-DIAL 632 
technique). This lidar with expected improvements will provide 633 
a unique data set to investigate the chemical and dynamical 634 
processes in the PBL and free troposphere. The spatiotemporal 635 
variance estimates derived from the o/one lidar observations 636 
will also be useful foi assessing the variance of tropospheric 637 
ozone captured by satellite retrievals 638 
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